Matter wave soliton collisions in the quasi one dimensional potential 
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We consider soliton solutions of a two-dimensional nonlinear system with the self-focusing nonlin- 
earity and a quasi- ID confining potential, taking harmonic potential as an example. We investigate 
a single soliton in detail and find criterion for possible collapse. This information is then used to 
investigate the dynamics of the two soliton collision. In this dynamics we identify three regimes 
according to the relation between nonlinear interaction and the excitation energy: elastic collision, 
excitation and collapse regime. We show that surprisingly accurate predictions can be obtained 
from variational analysis. 
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I. INTRODUCTION 

Solitons are universal entities in the nonlinear science 
and interactions between them are perhaps the most 
fascinating features. Understanding soliton collisions is 
both of fundamental interest and of importance to its 
wealth and proposed applications. If the system is in- 
tegrable, the collisions are elastic An example of 
such a system is ID soliton in attractive nonlinear Kerr 
medium. However, even if we weakly perturb nonlinear 
Schrodinger equation, we can observe a chaotic character 
of two-soliton collisions @, In the general case, it is 
necessary to resort to detailed numerical calculations for 
predictions. Stegeman and Segev [J] have introduced a 
general classification of collisions into two categories: co- 
herent and incoherent. Coherent interactions occur when 
the nonlinear medium can respond to interference effects 
that take place when the beams overlap. Incoherent in- 
teractions, on the other hand, occur when the relative 
phase between the soliton varies much faster than the 
response time of the medium [f|. In Ref. [1] and for 
example in Ref. Q one can find description of spatial 
soliton collisions in the Kerr medium. Collisions between 
spatiotemporal solitons of different dimensionality in a 
planar waveguide were considered in Here we also 
study collisions in the Kerr medium, which has infinitely 
short response time, but we are dealing with a pair of 
two dimensional solitons, both of the same dimensional- 
ity. We investigate head on collision of solitons moving in 
quasi ID confining potential. Our numerical results were 
obtained for the harmonic potential, but they generally 
apply to all quasi ID confining potentials that posses at 
least one bound state. We would like to point out that 
multidimensional solitons were extensively studied in the 
case of Bose Einstein condensates, including beautiful ex- 
periments @] . Some of the applications of solitons in the 
condensates was discussed in Ref. @. 

Within the framework of current publication the stat- 
ics and the dynamics of solitons in quasi ID potentials 
are described within the variational approximation and 



compared with full numerical simulations. We find for- 
mulas for the widths and chemical potential and discuss 
quasi ID limit and the dynamics of the collapse. Similar 
considerations were presented by Salasnich et al in refer- 
ence [l(| in the case of Bose-Einstein condensate [ll|, [H[ 
and earlier for the optical beam propagation by Li et al 
[l3| . Analogous study of the stability of gap solitons were 
presented in [l4| . Finally we address the problem of col- 
lisions, in various regimes, from elastic collision regime, 
through the domain where the transverse excitations oc- 
cur, up to the collapse during the collision. The main 
result of our study is the analysis of the collapse during 
the collision. We find it somehow surprising how good 
estimate can obtained from a simple variational model. 
We would like to point out that the deviation from one 
dimensionality in stationary properties and collisional 
dynamics of matter-wave solitons was recently investi- 
gated in [l5| using an effective one-dimensional Gross- 
Pitaevskii equation that includes an additional quintic 
sclf-focusing term. 2D soliton collisions were also studied 
experimentally in the case of spin wave envelope solitons 
in Ref. 0. 

II. QUASI-1D APPROXIMATION 

In this paper we consider two dimensional system with 
attractive nonlinearity, with an additional external quasi- 
1D potential. This system is described by 2D nonlinear 
Schrodinger equation (NLSE) 

iy t = ~(y xx +y yy )+vy-\2 D m 2 y, (i) 

where A2D > 0. We assume that the wavefunction is 
normalized to N. All the calculations were performed 
for the specific case quasi ID harmonic confining poten- 
tial V(x,y) = ujx 2 /2, see Fig. (JT|). Our predictions and 
physical properties derived here apply all the systems, 
two- and three dimensional, which have an external po- 
tential, confining in all but one dimension (in literature 
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FIG. 1: Quasi ID potential with harmonic confinement. 



they are sometimes called "potentials with transverse 
confinement") that support at least one bound state in 
the transverse direction. From the existing literature on 
that subject we would like to acknowledge three examples 
that are closely related to our study: l)Li et al [l3[ con- 
sidered the simple variational model of spatial solitons 
in planar waveguides, 2) The existence of 3D solitons in 
the inhomogeneous medium with harmonic potential was 
investigated by Raghavan [17], 3) The dynamics of sys- 
tems under transverse confinement were studied in BEC 
by Salasnich et al [Tpj . 

In the case of BEC, the equation (fTJ) can be derived 
from the full three-dimensional Gross-Pitaevskii equation 
for the wavefunction in physical coordinates ^(fjt) 



im f = V 2 ^ 

* 2m 



C/(f)§-A 3D |^| 2 *, 



(2) 



where X-^d = — 47ra s fi, /m, m is the atomic mass, wave- 
function norm N is the number of atoms, and a s is 
the scattering length. By assuming strong harmonic 
confinement lo z in the z direction, U(i) = muj 2 z 2 /2 + 
U±(x,y), that prevents excitation of higher modes of 
the trap in this direction, we arrive in dimensionlcss 
Eq. ([T]) after performing rescaling according to [x, y) — 
(x,y)/x , t — (Ti /mxpt, V = (mxl/fi 2 )U, X 2 d = 
(m/h 2 )y/ muj z /2nh\3D and \& = xq^>. Here xq is an 
arbitrary scaling parameter. 

In all the nonlinear systems with quasi ID potentials, 
as long as the transverse part of the potential supports at 
least one bound state, we expect, in certain range of the 
strength of the nonlinearity, to find 2D soliton solutions. 
More precisely, in such systems for small values of non- 
linearity we find stable solitons, but there is a threshold, 
a critical value of nonlinearity at which catastrophic self- 
focusing occurs. When we approach this critical value 
our soliton turns into the Townes soliton [l8[ , and above 
this value no stable solitons are available. 

As we mentioned above, to focus our attention in what 
follows we concentrate on the model with harmonic con- 
finement. In the limit when the trapping frequency be- 



comes large we expect that the energy associated with 
the nonlinear interaction becomes negligible in compari- 
son with fiLO. In this case our 2D soliton will become a 
product of a ground state in the direction of the confining 
potential and ID soliton in an unbound direction. The 
following calculation, based on the separation of vari- 
ables, confirms this statement. 

If we assume that our wavefunction ty(x,y,t) can be 
presented as a product of 



(3) 



where 4>(x) is normalized to unity, we can rewrite equa- 
tion dm 



(x)i>(y,t)} t = 



1 r 



~<f>xx(x) + ^U 2 X 2 (f>(x) 



~ip yy (y, t) + A 2D |^(y, t)\ 2 ^(y, t)\4>(x)\< 



cj>{x), (4) 



and if we assume that 4>(x) is a ground state of the trap- 
ping potential we obtain 



4>t - i^ 2 ^ 



<£= -^Vv> + a 2dM 2 W| 2 0. (5) 



Upon multiplying both sides of the equation by <p*(x), 
integrating over x and neglecting constant term we ob- 
tain one dimensional (NLSE) for ijj — t/j exp(— cu 2 t/2) (see 
Eq. (O below), with effective nonlinearity equal to 



Aid — A 



2D 



dx\c/)(x)\ 4 



(6) 



We will call this regime a quasi ID limit. Notic e, that i n 
the case of harmonic potential J dx\(j>(x)\ A = ^Jus/{2n). 

One dimensional NLSE is fully integrable. Conse- 
quently, we expect that two solitons, that are coming 
into collision will asymptotically reassume their original 
shape after collision. If the nonlinear energy becomes 
comparable with the excitation energy of the confining 
potential, we expect to see growing deviation from ID 
dynamics. It is the main purpose of this presentation to 
study this deviation. 



III. VARIATIONAL APPROXIMATION 



A. ID soliton 



The nonlinear Schrodinger equation in ID case reads 



ilpt = -\^yy - AidIV'IV, 



(7) 



where we assume that norm of the wavefunction is equal 
t° /_ IV'I 2 ^?/ = N . Within the variational approxima- 
tion [2(| , instead of solving Eq. ([7]) we introduce La- 
grange density 
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(8) 
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To mimic the dynamics we will introduce a trial function 
of the form (variational Ansatz) 



1>(y,t)=A(t)exp[-Y 



1 



v(ty 



- ib(t) 



+ #(*) , (9) 



with variational parameters: amplitude A(t), chirp 
b(t), overall phase <p(t) and width V(t). By substi- 
tuting Ansatz ^ into our Lagrange density and inte- 
grating over y we obtain the Lagrange function L = 



-4>v - 



bV 
4 
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iV 
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Our varia- 



tional approximation restricts the set of available solu- 
tions to the certain class of functions. It will lead to the 
Euler - Lagrange equations of the reduced Lagrangian 
shown above 
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with additional one, which we can interpret as a con- 
servation law (first integral) A 2 V = const, and which is 
related to the norm of the trial function 



\ip\ 2 dy = A 2 Vy/^=N 



A 



1 N 



(11) 



In this formalism solitons correspond to stationary solu- 
tions of Eq. (fit)]) i. e. we assume V = V = 0. These 

and 6 = 0. 



/2tt 
NXil 



conditions can be satisfied when V = 
Notice that we can evaluate soliton eigenvalue. In non- 
linear optics this eigenvalue corresponds to the soliton 
wavevector and in the theory of BEC it becomes chemi- 
cal potential. We write the solution in the form 



ip(y, t) = $(y) exp (-ifiiDt) ■ 



(12) 



We use 2D Gaussian Ansatz 



V(x,y,t) = A(t)exp 



x exp 1 



2W(t) 2 2V(t) 2 



^t) + -(b(t)x 2 +c(t)y 2 ) 



with variational parameters A(t), 4>(t), b(t), c(t), W(t) 
and V(t). In analogy with what we presented above we 
obtain (after integrating over x and y) the Lagrangian 
and look for the stationary solutions (W = W ~ and 
V = V = 0), which occur when the following conditions 
are satisfied: c = 0, 6 = and 



V = 



2tt 



\2DN 



:W, and W = 



At: 2 



Att 2 l 



(16) 



First significant observation is the clear evidence of col- 
lapse in our model. Note that equation for W can not be 
satisfied when X2DN > 2ir. Close to this critical point, 
both widths become equal and tend to zero. 

To obtain the value of the chemical potential we sub- 
stitute the width obtained in Eq. (|16p into equation for 
4> and get 



= 



2n 2 -X 2 D N 2 



uj \-2rj 2 



2tt ^Att 2 - X 2 2D N 2 2 ^1 



(17) 

where n = (\2dN)/(27t). In the analogy with ID case we 
can write ^(x,y,i) = e~ l/ - l2Dt $(x,y), Similar variational 
approach can be developed for the case of the quasi ID 
square well potential. The only difference in the func- 
tional form of the effective Lagrangian, and therefore 
also in the equations of motion, would be in potential 
term, which contains uj,. Finally we would like to point 
out that the analysis presented above can be used in the 
linear limit, when the waveguide mode structure can be 
predicted with satisfactory accuracy, see [T3 |. 



To evaluate we substitute V into the equation for 
phase, and obtain 



-fJ-iD 



5Xj D N 2 
8tt 



(13) 



B. 2D soliton 



In this case we start with 2D nonlinear Schrodinger 
equation with harmonic potential 

i¥ t = ~[* XX + ^yy) + ^Vtf - A 2D |*| 2 *, (14) 

with normalization/ \$>\ 2 dxdy — N. The 2D La- 
grange density is 



c = - [i(* t ** |ig 2 - 

-|vFJ 2 - w 2 x 2 |v]/| 2 + A 2 d|'I'| 4 1 



(15) 



C. Collapse of the 2D wavefunction 

We now consider the case when n > 1, i. e. in the 
regime where we expect the wavefunction to collapse. As 
we mentioned above, close to the collapse both widths of 
our solution, W and V become even. Hence, in the crude 
approximation, to describe the dynamics of the collapse 
we can assume axial symmetry. Upon neglecting har- 
monic potential contribution (which is negligible during 
the collapse) we obtain 



W 



1 



V 

w 3 



1 — 7] 

~W~' 



(18) 



If we solve the Eq. (|T8|) with initial conditions W(0) = 
Wo, W(0) = 0, which correspond to the dynamics origi- 
nated from some unstable state, we obtain 



(19) 
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FIG. 2: Comparison of the 2D soliton cross sections (dots) 
(in x and y planes) with quasi ID approximation (continuous 
curves), a) and b) corresponds to the frequency ui = 40. We 
see a perfect agreement between ID soliton and cross section 
of 2D soliton along y axis and ground state of the harmonic 
potential and 2D soliton cross section along x axis, c) and d) 
were obtained for smaller frequency ui = 0.5. Some deviation 
from quasi ID approximation can be observed. 



The main conclusion from this simple calculation is that 
the collapse of the wavefunction occurs within the finite 
time (except when r) = 1, on the threshold for the col- 
lapse). The collapse time is equal to 



(20) 



Collapse occurs on the timescale that is proportional to 
Wq , but what is more important it occurs the faster the 
the higher norm of the wavefunction is. This will be 
important in the next section when we discuss collapse 
during the soliton collision. 

To relate our result to experiments, we calculate 
the critical atom number necessary to observe collapse 
in BEC. We consider a 85 Rb condensate in a highly 
anisotropic trap configuration with oj z = 2?t x 350 Hz, 
Lu y = 2tt x 55 Hz, and a shallow confinement in the 
x direction. This configuration can be realized using 
an optical dipole trap 1211 . For the scattering length 
value of a s = — 15 ao [12j the critical number of atoms 
is N cx w 10 3 , and the typical soliton dimensions close to 
the collapse threshold are of the order of several /zm. 



IV. NUMERICAL RESULTS 

A. Quasi-ID limit 

We now consider the case when value of A2D (see 
Eq. ©) tends to zero, but the value of Am and the 




FIG. 3: Widths of the 2D soliton obtained from variational 
approximation, as a function of parameter A2D-/V (solid lines). 
Marked with dashed line are the corresponding values of ID 
soliton width and the width of the ground state of the har- 
monic ID potential (horizontal line). 
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FIG. 4: Widths of the 2D soliton obtained from direct numer- 
ical simulations, as a function of nonlinear parameter \2dN . 



norm remain constant. This corresponds to the limit 
lu — ► 00. By substituting 2D nonlinear coefficient A2D = 
Aid \/2it/u) into the expression for 2D eigenvalue [i2D in 
Eq. ([XT]), we obtain 



uj 3A 2 1D N 2 



8tt 



2 +A*ir> ■ 



(21) 



This way we show that 2D eigenvalue, in the limit of high 
frequency, consists of ID eigenvalue and the ground state 
energy of harmonic oscillator. It is a very intuitive result, 
since in this limit 2D soliton becomes a product state of 
ID soliton and the ground state of the harmonic potential 
in the transverse direction. Now we examine the widths 
of 2D soliton presented in Eq. (JTHJ) . If we assume to — > 00 
one of the widths W ~ and V -> V2tt/(A 1d N). 

This is exactly the value we obtained for the ID soli- 
ton. In conclusion, in the limit considered here, both 
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FIG. 5: Comparison of the chemical potential /j,2d obtained 
using quasi ID approximation (Eq. (|21[) (dashed line) and 
full 2D variational approximation (Eq. (I17|l (solid line) as a 
function of nonlinear parameter \2dN. Notice that the region 
between the collapse and the point up to which quasi ID 
approximation holds is very narrow. 



chemical potential and the width obtained from varia- 
tional approximation take proper quasi ID values. In 
Fig. [5] we present both cross-sections of 2D soliton in x 
and y planes. One of the cross sections is compared with 
ID soliton and the other with the ground state of the 
harmonic potential. The comparison is made for two dif- 
ferent values of frequency to. For larger value of u we ob- 
serve perfect matching and when u> becomes smaller, we 
can see slight deviations from quasi ID approximation. 
In Fig. [3] we present 2D soliton width obtained from vari- 
ational approximation, as a function of parameter A2D-/V. 
For comparison we also included a dashed line repre- 
senting a width of the ID soliton, and straight dashes 
horizontal line corresponding to the width of the ground 
state of the trapping potential (Gaussian). In Fig. 2] the 
same comparison is made for the case of numerical so- 
lutions. In this case we defined the widths numerically, 
as (VjWO = -^((M): (|y|)); i- e - the the mean value 
of the modulus of the coordinate in this direction. We 
clearly see that for small \2dN 2D soliton width well ap- 
proximates the equivalent ID value, and when \2dN ap- 
proaches critical value, both widths become equal, while 
tending to zero. In this limit, the energy of the trap- 
ping potential is much smaller than kinetic and nonlinear 
energies, hence our solution becomes practically identi- 
cal as that of Townes soliton [HI]. Finally, in Fig. [5] we 
plot a chemical potential as a function of \2dN obtained 
within quasi ID approximation (Eq. (|2ip ) and that ob- 
tained from full 2D variational analysis, (formula (|17|1 ). 
Notice that when the \i2D is of order of the excitation en- 
ergy in the harmonic potential quasi ID approximation 
breaks down. Also on this figure we can see a clear indi- 
cation of the collapse region, when A2D -/V approaches the 
value of 2ir, and the fact that region between the collapse 
and the point up to which quasi ID approximation holds 




FIG. 6: Boundaries of the collisional collapse region for the 
trap with the frequency w = 4, obtained from variational 
approximation (with Gaussian - upper continuous curve and 
hyperbolic secant trial functions - lower continuous curve) and 
obtained from direct numerical simulations (crosses) . Solitons 
colliding with velocities below those marked with crosses will 
experience collapse during their interaction. 



is very narrow. 



B. Soliton collisions 

Now we investigate collision of solitons that are mov- 
ing along the quasi ID harmonic potential presented in 
Fig. [U Soliton collisions in the one-dimensional NLSE 
were studied in many contexts, including nonintegrable 
dynamics of vector solitons (2^]. In the case of BEC, the 
quasi- ID collisions were investigated in Rcf. in the 
framework of the nonpolynomial Schrodinger equation 
(NPSE). The results obtained here can be easily general- 
ized for the quasi ID square well potential. The solitons 
that we used in the simulations were identical, and had 
equal and opposite velocities. 

We start with a pair of solitons, each of which is sepa- 
rately a solution of Eq. |T]). The initial width along free 
direction we denote by Vo and the in the transverse di- 
rection by Wo . We assume that the norm of each soliton 
as equal to N/2 and they travel with velocity v. The 
nonlinear interaction is described by the parameter \2D- 
We will introduce here a simple estimate of collapse dur- 
ing the collision based on variational approximation. It 
is well known that estimates based on variational analy- 
sis depends on the choice of variational basis, up to the 
multiplicative factor [23|. Hence, to improve the qual- 
ity of our predictions we will consider the estimates ob- 
tained using Gaussian functions and hyperbolic secant 
functions. Lets first concentrate on Gaussians. The first 
parameter that describes collision is soliton interaction 
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time. This is merely the time when both wavefunctions 
have a significant overlap. If we assume that significant 
interaction appear as long as wavepackets are separated 
not more than FWHM, we obtain interaction time equal 
to 



2jV a _ AjWq 
v rjv 



(22) 



where 7 ~ 1.178. Notice that the width ratio differs 
from the one given in Eq. (|16p . because soliton norm is 
now N/2 instead of N. Here we also assume that during 
this time interaction is almost constant and is taken as 
that of full overlap (both function on top of each other) . 
Now we turn to the condition for the collapse. As we see 
from Sec. IIII CI collapse is expected to occur for 77 > 1. 
This condition, in the case of collision, should refer to 
the situation when wavepackets fully overlap. On the 
other hand soliton with the norm N/2 can only exist if 
77 < 2. Taking both conditions into account we conclude 
that the region of interest is 1 < 77 < 2. We estimate a 
time of collapse for solitons in the above interval to be 



W 2 



(23) 



Notice that here Wo is the initial width of each of the 
colliding partners and given by 



W a = 



4 - i] 2 
Aoj 2 



(24) 



From the above discussion it follows that there should be 
a critical velocity, above which during the collision the 
collapse will take place. This condition can be obtained 
by comparing time of collision and time of collapse 



W 2 



r]v c y/rj - 1 
which after some algebra will give 



(25) 



v c (v) 



47^4 (ry-1) 2 



n 



(26) 



If the velocity of solitons exceeds this critical value, there 
is not enough time during the collision to complete the 
collapse. If the velocity is lower than this value, during 
the overlap time collapse can fully develop. Hence the 
function v c (r]) marks the boundary between regions of 
collapse and no collapse. 

As we mentioned above, position of this line depends 
on the choice of the trial functions, which in the case just 
described were Gaussians. Had we chosen hyperbolic se- 
cant functions instead of Gaussians, we would have ob- 
tained slightly different result. First, we find that in this 
case the collapse for single soliton appear for the value of 
nonlinearity (A2D-/V) equal to 6 instead of 2tt. Hence, the 



parameter 77 should be defined as (A2_d-/V)/6, and after al- 
gebra, very similarly to what we presented for Gaussians, 
we obtain the condition for the critical velocity as 



v c (v) 



87^ 4/ (v — I) 2 

V^v V 
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(27) 



In Fig. [6] we show curves representing v c as a function 
of \2dN, obtained from Eqs. (|2"6")) and (|2"7|) . indicated 
with solid lines, together with the results obtained from 
direct numerical simulations (crosses). The agreement 
between the variational predictions and direct numerical 
simulations is very good. It is somehow surprising. It 
was shown in Reference [lij that the variational anal- 
ysis of the collapse dynamics does not work very well 
close to the threshold. Nevertheless it seems that for our 
crude estimate it is sufficient. All the crosses are falling 
between curves obtained with Gaussian and hyperbolic 
secant trial functions. Solitons colliding with velocities 
below those marked with crosses will experience collapse 
during their interaction. In conclusion, we see that varia- 
tional analysis predicts correctly and accurately the onset 
of catastrophic collapse during soliton collision. 

The full picture of the soliton collision that emerges 
from our numerical studies is the following. Imagine we 
keep the norm of each soliton and their velocity constant. 
Then, depending on the strength of the confining poten- 
tial (u> in the case of harmonic potential and the width 
in the case of quasi ID square well potential) we can 
clearly identify three regimes according to the relation 
between nonlinear interaction and the excitation energy. 
In the first regime, when the distance between ground 
state and first excited state of the confining potential is 
much bigger than the interaction energy, system is prac- 
tically integrable, hence solitons pass though each other 
and restore their original shape after the collision. When 
both energy scales become comparable we observe exci- 
tations in the transverse directions, corresponding to the 
transition between ground state and first coupled excited 
state. Finally, we observe catastrophic collapse during 
the collision. Behavior of solitons in these three regions 
is illustrated in Fig. [7] Picture a) is characteristic for 
the first region described above. It was obtained for high 
frequency u> — 20 and small nonlinearity A2D = 2. Soli- 
tons appear to be practically one dimensional, and they 
asymptotically (after the collision) restore their original 
shape. In the case b) lu = 2 and nonlinearity is the 
same as in previous case and nonlinear interaction be- 
comes comparable with excitation energy. Some of the 
population will be transferred to the excited state and we 
observe the beats of the frequency 2u> in the transverse 
direction. It can be detected for example by calculating 
mean square radius in the transverse direction, shown in 
the inset on the right hand side. Finally, in picture c) 
where A215 = 3 and u> = 2, total population is above 
the critical value and collision time is sufficiently long - 
we observe collapse during the collision. Analogous ef- 
fect was previously observed in optics, when two spatial 
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FIG. 7: Gallery of collisions. From quasi ID to collapse. Soli- 
ton velocity v — 2.5. Panel a) high frequency (u — 20) and 
small nonlinearity X2D = 2 - solitons appear to be practically 
one dimensional, and they asymptotically restore their origi- 
nal shape. Panel b) small frequency u> — 2 and nonlinearity 
the same as in previous case - nonlinear interaction is com- 
parable with excitation energy. We observe oscillations with 
frequency 2cj (inset on the right hand side shows mean square 
radius in the transverse direction as a function of time) . Panel 
c) frequency the same as in case b), but higher nonlinearity 
A2D = 3 - we observe collapse during the collision. 



solitons collided in planar waveguides, see for instance 



V. CONCLUSIONS 

In conclusion, we analyzed the stability and collisions 
of quasi ID solitons in the confining potential, both 
numerically and within variational approximation. We 
showed that variational approximation constitutes an ex- 
cellent basis for estimating physical parameters of quasi 
ID solitons. The advantage of using this approximation 
is that one obtains analytical formulas, and can predict 
soliton behavior in different limits and under various cir- 
cumstances. One of the examples is the threshold of 
the collapse and its dynamics. Comparing time of col- 
lapse and collision time we could predict the result of 
the soliton collision and identify three different regimes 
according to the relation between nonlinear interaction 
and the excitation energy. In the first regime, when the 
distance between ground state and first excited state of 
the confining potential is much bigger than the nonlin- 
ear interaction, system is practically integrable, hence 
solitons pass through each other and restore their orig- 
inal shape after the collision. When both energy scales 
become comparable we observe excitations in the trans- 
verse directions, corresponding to the transition between 
ground state and first coupled excited state. Finally we 
observe catastrophic collapse during the collision. In the 
following paper we will present similar considerations for 
the case of Gap solitons. 
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